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On the determination of the decay constant of a radio- 
active substance 


By NgAu Hoe 


In this paper only the case will be treated when the decay of a radioactive 
substance is recorded by means of a counting device. 

The method most frequently applied for determination of the decay constant 
(or what is equivalent to this, the half-life) is a graphical analysis of the decay 
curve. Numerical methods have, however, also been proposed and used. In 
this note we are going to derive some formulae of this kind by means of statistical 
methods. We have chosen a starting point which, for the case treated here, 
leads to relatively simple calculations. The reader, who is interested in the 
formulae derived, but not in the derivations, may consult the summary. 

We shall confine ourselves to the case of a radioactive source consisting of 

one single radioactive substance. In the calculations is further assumed that 
no background radiation is present. This is in most cases to be considered as 
a first approximation. The results obtained when zero background is assumed, 
are, however, quite applicable also in the case when the background is small 
compared to the intensity of the radiation from the source. 
Consider one of the radioactive nuclei in the source, which at some arbi- 
trarily chosen zero point on the time axis has not yet disintegrated. The 
probability that it will disintegrate between ¢ and ¢ + dt is then, according to 
the well known probability law governing radioactive decay 


P(t)dt = ke dt. (1) 


The parameter k of this formula is the decay constant of the radioactive 
substance. It is connected with the half-life 7, by the equation 


kT, =1n2. 


Our experiment consists in continuous observation of the radiation emitted 
from the radioactive source during some time interval, say from 0 to 7. For 
the sake’ of simplicity we shall first consider the case when all disintegrations 
taking place in the source in this time interval are registered, a situation which 
does not occur in practice. 

As we observe only those disintegrations which take place between 0 and 7, 
we may choose as a starting point for our calculations, not the distribution 
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with the probability density P(t) of equation (1), but the conditional distribu- 
tion of t relative to the hypothesis 5 


Mee bcs, 


We then get a truncated distribution with the probability density 


p(t) = ape () 


in the above mentioned interval for ¢. This distribution is a one-dimensional — 


one with one statistical parameter k which has to be estimated from the ex- 
perimental data. 

Before starting the derivation of formulae for estimates, were are going to 
modify our starting point so that it is more in accordance with real experiments. 
Equation (2) is derived on the assumption that all disintegrations occurring in 
the source are observed. In almost every experiment performed in practice 
this is not the case, partly because only a fraction of the radiation emitted 
from the source hits the counter, and partly because the efficiency of the 
counter is less than 1. 

It can, however, easily be shown that equation (2) is valid also in the more 
general case when only a fraction of the disintegrations is registered, and that 
this equation can accordingly be used as the basis of the calculations also in 
this case. 

Consider a statistical experiment a with a given probability density p (é), 
Then consider a composite statistical experiment, where first the experiment a 
is performed and then a new experiment which decides whether the first ex- 
periment is to be counted or not. The probability that it is to be counted is 
some given number. Then the results obtained by this composite experiment 
are also distributed with the probability density » (t). 

From these considerations it follows immediately that distribution (2) is also 
valid for the counts registered in a counter which is exposed to the radiation 
from a decaying radioactive substance. This is of course true only if the loss 
in the counting device on account of its finite resolving power is negligible. 
We shall base our calculations on the assumption that this loss is negligible. 
Corrections may then, if necessary, be made in the formulae derived. 

We now consider our physical measurement which consists in the observation 
of all counts of the counter in the time interval from 0 to 7. Let the number 
of counts be denoted by m. The measurement may then be considered as m 
independent observations of a statistical variable ¢ having the probability 
density p(t) as given by equation (2). Our problem is then to find an estimate 
for the parameter k& from these observations. For this purpose we are going 
to use the method of maximum likelihood. 

Suppose that the counts have occurred at t, to t, + dt, ty to ty + dt, ..., tm to 
tm + dt. The probability for this result is 


k m =e ty 
é 


S, (dt)” =m! (, oT =1 (deym, (3) 
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According to the method of maximum likelihood we have to choose as our 
estimate for & the value of k, which makes S, as great as possible. This value 
is In this case found from the equation 


In this equation we replace k by x, in order to make clear that we are 
dealing with an estimate for k and not the parameter itself. We write x, in 
order to distinguish between this estimate for & and another one x, which is 
derived later. 

For the determination of x, we find the following equation 


- m 
1 1 fi 


Tu, eF—] mT 


(4) 


This formula was first found by Prrerts [1] who used the method of inverse 
probability, and it was later derived by Hote [2][3] who used the method of 
maximum likelihood. 

It has been shown in a preliminary communication [8] that it is possible to 
find the exact sampling distribution of x,. Here we shall only give the result 
of the calculation. 

We first consider the statistical variable 


and denote its probability density by gm(u). The problem to find the distribu- 
tion of u, is closely connected with the problem of finding the distribution of 
the sum of m observations of a statistical variable with a rectangular distri- 
bution. 

For gm(u) we find 


BT*\* 
Jm (u) = (, | 8d ed OD) (5) 
with 
= Jie m=—t —— — ]\m—1 (ss 9\m—1 __ 
The summation is continued as long as the arguments u,uw—1,w—2... are 


positive. The distribution of x, may then be found by means of the equation 


ses 1g 
Coe ONT elm —] 


The use of the exact distribution is, however, complicated, and as a rule one 
is mainly interested in some simple measure of the dispersion of the estimate. 
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The estimate x, has been found by means of the method of maximum like- 
lihood, and it is known that under certain general conditions (see for instance [4]) 
an estimate found by this method is an asymptotically efficient estimate. This 
means in our case that the estimate x, may for practical purposes be con- 
sidered as an efficient estimate. We can therefore apply the formula for the 
mean square deviation of an efficient estimate to find the mean square devia- 
tion o7 (x1) of 2%, (see for instance [5]) 


(1) = — aay () 


0? In S, 
os (— e 


By means of equations (3) and (7) one easily finds 
e 


o (x) ag i 
oe (1 zs ay ®) 
m (et — Ne 


So far we have assumed that the time for every count is recorded, and in 
the determination of the decay constant use has been made of all this know- 
ledge about the counts. This method is, however, as a rule, not applied in 
practice. 

If the source is observed continuously from 0 to J’, it is usual to divide 
this time interval into 7 intervals of equal length 


and to observe the number of counts in each of these intervals. This corre- 
sponds to a grouped set of observations of the statistical variable ¢ of equa- 
tion (2). 

We denote the number of counts in these intervals by x, %,...a,. From 
these data we are then going to find an estimate for the decay constant. 

This problem can be given the following form: Consider the number of a 
time interval as a random variable of the discontinuous type, capable of as- 
suming the values 1, 2,..., 7. The probability that this variable takes the value 
2 is p;. One then has 


pi= | p(t)dt 
(i-1)1 
The integration gives 
Lee. 
Di = Fane eee: (9) 
The probability that our observations give the result 21, v2, ... a, is accordingly 
m! , 
ies 2 py & av, 
Se £1! 2%! enn a bi ae te 
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where 
M = Ly + yg tes? + Bp. 


If (9) is inserted for p,, pg, etc. one obtains 


! — e-kl\m 
Ss ade = ( l~e em FU ae tay +++ + (1p), (10) 
ns 


ary rae 1—e*t 


The application of the method of maximum likelihood then leads to the same 
procedure as in the case first treated, and we get the following equation for 
the determination of the estimate xg: 


> 1 Yo ee + ie, ee (pe 1a 
ga —] eT —] m A 


_ UF Aa, ti The ibaa ah 
m 
This equation was first found by Hox [2]. If the loss of impulses in the 
counting device is not negligible, one should insert for x1, 2% etc. the numbers 
obtained after correction for this loss has been made. | 

For the practical application of formula (11) the function 


ih ear 
e&— ] ei — | 


F, (q) 


has been tabulated for some values of r and gq, and the results are given in 
Table 1. 

It is also possible in this case to find the exact sampling distribution of 
the estimate. This problem is solved if the sampling distribution of 


Uf =, bag a 2° He Tay 


can be found. This problem is, as has been shown [6], closely connected with the 
well known classical problem of De Moivre. We are not going to give details 
of the derivation of the distribution of y here. If the probability for the value 
y be denoted by x(m,y), it may be shown that 


a(m, y) = A™ a’ G(m, y) (12) 
where -, 
ee’ — 1 
NS 1—e-*T 
a = etl 


Lely ed \e Yom" yr m ae 
Tod bei end Ma eG lec OD 
When y Sr +™, the formula has only one term. When 2r +m2y>r+™m, 


it has two terms etc. 
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q Ps Fo ys F16 Fx9 
0.125 2.8546 3.4950 4.0638 5.0062 5.7219 
0.130 2.8297 3.4568 4.0106 4.9189 5.5985 
0.135 2.8049 3.4190 3.9580 4.8329 5.4777 
0.140 2.7802 3.3814 3.9057 4.7482 5.3596 
0.145 2.7556 3.3441 3.8540 4.6649 | 5.2441 
0.150 2.7311 3.3070 3.8028 4.5829 5.1313 
0.155 2.7067 3.2702 3.7520 4.5022 5.0210 
0.160 2.6824 3.2336 3.7018 4.4228 4.9134 
0.165 2.6583 3.1974 3.6521 4.3448 4.8084 
0.170 2.6343 3.1614 3.6030 4.2682 4.7060 
0.175 2.6104 3.1256 3.5543 4,1929 4.6061 
0.180 2.5866 3.0902 3.5063 4.1190 4.5087 
0.185 2.5630 3.0551 3.4587 4.0464 4.4138 
0.190 2.5395 3.0202 3.4117 3.9752 4.3213 
0.195 2.5161 2.9857 3.3653 3.9053 4.2312 
0.200 2.4929 2.9515 3.3194 3.8367 4.1435 } 
0.205 2.4698 2.9176 3.2741 3.7695 4.0581 
0.210 2.4469 2.8839 3.2294 3.7036 3.9750 
0.215 2.4241 2.8506 3.1852 3.6390 3.8940 
0.220 2.4014 2.8177 3.1416 3.5757 3.8152 
0.225 2.3789 2.7850 3.0986 3.5137 3.7385 
0.230 2.3566 2.7527 3.0562 3.4530 3.6639 
0.235 2.3344 2.7207 3.0143 3.3935 3.5913 
0.240 2.3123 2.6890 2.9730 3.3352 3.5207 
0.245 2.2904 2.6576 2.9322 3.2782 3.4520 | 
0.250 2.2687 2.6266 2.8921 3.2223 3.3852 
0.255 2.2471 2.5958 2.8525 3.1676 3.3201 
0.260 2.2257 2.5655 2.8134 3.1141 3.2569 
0.265 2.2044 2.5354 2.7750 3.0618 3.1953 
0.270 2.1833 2.5057 2.7371 3.0105 3.1354 
0.275 2.1624 2.4763 2.6997 2.9604 3.0772 


By means of the equation 


=m\1 + : aes. 
4 gh] * fa] 


one then gets the exact sampling distribution of xg. 

The use of the exact sampling distribution is complicated in this case too, 
and we therefore ‘prefer to use the expression found for the mean square devia- 
tion of x, when we assume that m is so great that x, may with sufficient ac- 
curacy be considered as an efficient estimate. This is justified for all practical 
purposes. | 

We now use formula (7) again after having written x, instead of x, and S» 
instead of S,. After some simplifications we find [7]: | 
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GO" (x2) 1 
2 2 kl 2 
k m( (kl)? e (kT) _ 


(et! — i (e*?—1)3 


(14) 


It is easily shown that formula (14) becomes formula (8) in the limiting case 
when r->ocoo: 10. 
For the practical use for formula (14) the function 


D; (q) = (@ Ie = Wy eal” 


has been tabulated for some values of ry and q, and the result is given in 


Table 2. 


: Table 2 

: 

q Dg | Dy D4 Di4¢6 Dy 

: 

0.125 12.8 8.4 6.0 3.64 Pail 

0.150 9.1 6.0 4.4 Qallo) 1.99 
0.175 6.9 4.6 3.4 2.19 1.66 

: 0.200 5.4 ail al 1.84 1.44 

: 0.225 4.4 3.0 DES) 1.61 L3l 

0.250 esait 2.6 220) 14s) iL Pall 

| | 0.275 332 Ba 3 1.8 LB 1S 


A value found for the mean square deviation in this way must be considered 
as a lower limit for this quantity, as there besides the purely statistical dis- 
persion of the x-values, will also occur fluctuations due to imperfections of 
the experimental arrangement which is used. 

Further it must be noticed that in most experiments a background radia- 
tion is also present. In such a case the method of maximum likelihood does 
not give formula (11) for x ,, but a much more complicated expression [2] 
which is not suited for practical application. If the background radiation is 
not too great, the estimate given by formula (11) gives a good estimate for 
the decay constant. For x; one must in this case insert the recorded number 
of counts minus the average number of background counts, which is as a rule 
known. 

In this case the mean square deviation given by formula (14) is only a 
lower limit even for an experimental arrangement functioning without error. 
If the background is low as compared to the intensity of the radioactive 
source at the end of the measurement, one may use (14) as an approximation 
to o7 (x2) and add a little to this value to correct for the approximation. This 
must also be done because of other errors not involved in the derivation of 
(14). This uncertainty in the estimation of o7 (x2) is, however, not so serious 
because if great precision is wanted, many activations are made. The final 
value accepted for x, is then the mean of a series of values, and o7(x,) may 
be estimated also from the dispersion of these values. The value found in this 
way together with the one based on equation (14) should then give a quite 
reliable result. 
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SUMMARY 


In this paper formulae for the numerical treatment of the measurement of 
the decay of a radioactive substance are derived. First the case is treated, 
when the time for all counts between 0 and 7 is registered. An estimate 
for the decay constant for this case is given in formula (4). The exact sampling 
distribution of this estimate is then given, and finally equation (8) gives an 
expression for the mean square deviation of the estimate. 

Then the case is treated which is of more interest to the experimentalist. 
In this case the decay of the radioactive substance is also followed continuously 
from 0 to TZ. This interval is divided into 7 smaller intervals of equal length 


= * and the number of counts in each of these intervals is recorded. An 


estimate 2 for the decay constant in this case is given by formula (11). The 
exact sampling distribution of the estimate is derived, but turns out to be too 
complicated for practical use. Finally the mean square deviation of the estimate | 
%, 18 given by formula (14). | 
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